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We study the LHC search constraints on models that extend the Standard Model (SM) with an 
inert, complex scalar electroweak multiplet, E, with isospin T = 5/2 (sextet) or T = 7/2 (octet) 
and identical hypercharge to the SM Higgs doublet. Imposing a global Z 2 symmetry under which 
E ^ — E, the lightest member of E is stable and we require that it be neutral (C^’’^) to avoid 
cosmological constraints from charged relics. Pair production of scalars by electroweak interactions 
followed by cascade decays to through W and Z emission produces signatures similar to those 
of supersymmetric electroweak gauginos, and we constrain the models by recasting a collection of 
such searches made with data from the 8 TeV run of the LHC. We find that there is no sensitivity 
from these searches to the compressed spectrum regime, in which the mass splittings between the 
lightest and heaviest states in E are less than about 20 GeV. In the remaining parameter space, we 
find significant exclusions for m^o,r ~ 80 — 180 GeV in the sextet model, and ?Ti^o,r ~ 80 — 120 GeV 
in the octet model. 
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I. INTRODUCTION 


The discovery of a Standard Model (SM)-like Higgs boson at the CERN Large Hadron Collider (LHC) [TJ[2] provides 
proof of the existence of weakly coupled scalar particles in nature. Extensions of the scalar sector of the SM are thus 
of great interest, particularly when they could be connected to the solution of other mysteries of the SM, such as the 
nature of dark matter (DM) in the Universe [3]. 

Extensions of the SM Higgs sector involving an additional “inert” scalar multiplet, the lightest state of which is 
stable and hence a possible dark matter candidate, have been well-studied in the singlet [4], doublet [5], and triplet [6] 
cases (for recent summaries of the experimental status of these models, see, e.g.. Refs. [7], [H] and ED- More recently, 
multiplets from larger representations of SU{2)l have been investigated in the context of dark matter [T0HI2]- 

In this paper, we expand upon the study in Ref. [12], which focuses on models in which the SM is extended by a 
scalar electroweak multiplet with half-odd-integer weak isospin, T = (n — l)/2, where n is an even integer that counts 
the number of complex fields in the multiplet. This scalar multiplet has weak hypercharg^U = 1 (the same as that of 
the SM Higgs doublet), and transforms as odd under an imposed global Z 2 symmetry, allowing for the lightest member 
to be a stable dark matter candidate. In particular, we focus on the models with weak isospin T = 5/2 (sextet, n = 6) 
and T = 7/2 (octet, n = 8)0 Larger complex scalar multiplets are disallowed by perturbative unitarity considerations 
involving the scattering of two scalars into two gauge bosons mi. Smaller multiplets, such as the T = 3/2 case, have 
previously been studied in the literature m- 

In these models, the large number of new scalar particles leads to a plethora of electroweak pair-production processes 
accessible at the LHC via an s-channel photon, lU, or Z. Because of the Z 2 symmetry, these scalars can decay only 
to lighter scalars within the large multiplet via the emission of VE or Z bosons, terminating at the lightest state 
of the large multiplet, which is stable and escapes the detector. These production and decay modes produce 
similar phenomenology to supersymmetric electroweak gauginos, and the most promising searches to constrain these 
models are similarly those with a combination of missing transverse energy {$t) leptons. We avoid signatures 
with too many high-energy jets in order to avoid large QCD backgrounds, such as ti. In particular, we have recast a 
number of existing searches for supersymmetric electroweak gauginos by both the ATLAS pRll25] and CMS [26ll28] 
collaborations which use data taken at a centre-of-mass energy of 8 TeV (LHC8). 

This paper is organized as follows. In Sec. [H] we introduce the two large-multiplet models first studied in Ref. m- 
We re-express the parameter freedom in the models in terms of physical masses and mass splittings, which are most 
physically relevant for the kinematics of pair production and decays. We also translate the theoretical and indirect 
experimental constraints on the models previously studied in Ref. [12] into this new parameter space. In Sec. |HI| we 
discuss the pattern of decay branching fractions of each of the scalars. In Sec. IV we recast a collection of LHC8 
searches for supersymmetric particles to constrain our models. We conclude in Sec. [V| Some mathematical details 
and simulation issues are relegated to the appendices. 


II. DESCRIPTION OF THE MODELS 


We consider two models that extend the SM through the addition of a single, large electroweak multiplet of complex 
scalars, E, which has quantum numbers under SU{3)c 0 SU{2 )l U{1)y of (l,n, 1), where n = 6 or 8 is the size of 
the multiplet. In these models, the most general gauge-invariant scalar potential that preserves a Z 2 symmetry under 
which E ^ — E is given by 


y(i>, S) = + Ai + A3$V“i> 


+ 


A4 I]^T“S + h.( 


+ o(s"), 


( 1 ) 


where 4> is the SM SU{2)l doublet. Here ^ = C4>* and E = CE* are the Higgs doublet and the large scalar multiplet 
in the conjugate representation, respectively. The conjugation matrix, C, is an antisymmetric n x n matrix equal to 
for the SU(2)i;, doublet, and is given in Appendix [A| for the n = 6 and n = S representations. The and 
matrices are the generators of SU (2) l in the doublet and n-plet representations, respectively. 

The term E^T^E can only be non-zero when n is an even number, i.e., for half-odd-integer values for the total 
isospin, T, of the scalar multiplet. Together with the constraint T < 7/2 (n < 8) for complex scalar multiplets that 


^ We normalize Y such that Q = Y/2. 

^ A recent study m of the renormalization group running of the couplings in such models found that the sextet model’s couplings will 
diverge at a scale below 5000 TeV, if the new scalars are at the 100 GeV scale. The octet model was not studied explicitly in Ref. m, 
but the stronger running due to the larger electroweak representation means that its couplings will diverge at an even lower scale. This 
implies either that the new scalars must be composite particles or that additional fermions must enter in such a way as to sufficiently 
modify the running of the scalar couplings of the large multiplet. 
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arises from perturbative unitarity of scattering amplitudes involving the large multiplet’s weak charge m, this limits 
the models of interest to the cases T = 5/2 (n = 6 ) and T = 1 j2 (n = 8)j^ For these two cases, the large multiplet is 
given in the electroweak basis by 

S(n=6) = (c+t c+t c+t c°, rt c-^)^, 

S(„=8) = (c+t c+t c+t c+t c°, c“t rt c^f ■ ( 2 ) 

Note that the conjugate of the charged state is written as which is not the same as C~^- 

When the A 4 term in Eq. § vanishes, the Lagrangian preserves a U{1) symmetry (models with such a U{1)- 
symmetric potential have been studied in Ref. m ). The A 4 term breaks the U{1) down to a Z 2 symmetry [ 12 ], and 
splits the complex neutral component of E into its real and imaginary parts, Re and Im 

with different masses. Furthermore, the A 4 term induces and controls the amount of mixing between the charged 
states with the same electric charge, and C~^*- The mass eigenstates are defined for Q > 0 in terms of a mixing 
angle aq such that 


= cos aq + sin aq ( , 

= - sin aq + cos aq , (3) 


with m^Q < m^Q. Since there is only one state with \Q\ = n/ 2 , it remains unmixed. Expressions for the mixing 

angles are given in Eqs. (7) and (14) of Ref. [ 12 ], and the details are summarized in Appendix [a| for completeness]^ 
For these models to contain a dark matter candidate, we require that the lightest (stable) member of the large 
multiplet be electrically neutral. This occurs only when jAsj < 2 |A 4 |. In addition, and without loss of generality, we 
choose the real part of to be the lightest member of the large multiplet; this constrains the sign of A 4 such that 
A 4 < 0 for the sextet model and A 4 > 0 for the octet. 

The masses of the physical states in terms of the mass of the neutral real particle, and the Lagrangian 

parameters A 3 and A 4 are given by [ 12 ] 


m^o,i = m^o 




+ ^4, 


to^+q = m?o,r 

ill 2 ^ 




n 


i(-l)t A 4 T ^g2A2 + („2_4g2)A2 


C +2 


[A3 + 2(-l)t+iA4] , 


(4) 


where the notation is such that the sign in m?Q forces the relation m^Q <m„Q. The physical scalars arising from 

ill 2 ^2 

the large multiplet always occur in the same hierarchy, given from lightest to heaviest by: 




/-i3 /■O 4 

^±± ^±3 /-±4 ^±3 ^ihih Tj^ /'O 4 
-^1 5-^1 5S) 9-^2 9-^2 9-^2 9 S 


(5) 


Sample mass spectra are shown in Fig. where the two left plots are for n = 6 and the two right plots are for n = 8 . 
For each of the two plots for each multiplet, we have fixed one of A 3 or A 4 and varied the other to illustrate the effect 
of varying these parameters. 

As seen in Eq. 0 - the states of E interact with the SM through a coupling to the Higgs doublet, as well as through 
their gauge-kinetic terms. The couplings to the photon and the Higgs boson are diagonal in the mass basis, and thus 
do not induce decays. However, off-diagonal vertices involving two different scalar mass eigenstates and a Z or VF 
boson do exist. Thus, the decays of the mass eigenstates in E to lighter members of E occur only through emission of 
a kF or Z, which may be off-shell depending on the mass splitting involved. This leads to a distinctive decay pattern, 
which is shown for the n = 6 and n = S models in Fig. |^for a typical parameter point. 

Reference m studied theoretical and indirect experimental constraints on the model parameters (m^o,r, A 2 , A 3 , A 4 ), 
which we apply here: 

(i.) The zeroth partial wave 2^2 scattering amplitude involving pairs of scalars must satisfy perturbative unitarity. 


^ The model with T = 3/2 has been studied in Ref. | 15| . 

^ We correct a typographical error in Eq. (A7) of Ref. for the generic expression for the mixing angle aq. The expressions in Eqs. (7) 
and (14) of Ref. |12| are correct. 
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FIG. 1: Sample mass spectra for fixed values of A 3 or A 4 . The left two plots are for the n = 6 model and show, from 
top to bottom, (dashed black), H 2 (dashed red), (dashed blue), (solid green), (solid blue), 

(solid red), and (solid black). The right two plots are for the n = S model and show, from top to bottom, 
(dashed black), (dashed red), (dashed blue), (dashed green), (solid brown), (solid green), 

(solid blue), (solid red), and (solid black). 



FIG. 2 : The decay patterns of the n = 6 (left) and n = 8 (right) models for a typical parameter point. The upper 
branch of states includes the and while the lower branch includes and ^ each sorted by 

electric charge. Decays are shown as diagonal (IF emission) or vertical {Z emission) arrows. On-shell decays are 

indicated by a solid line and off-shell decays by a dashed line. 


|Reao| < 1 / 2 . This results in constraints for the sextet model of IA 2 I < 6.59, jAsj < 8.48, and IA 4 I < 4.25, and 
for the octet model of IA 2 I < 3.10, IA 3 I < 5.46, and IA 4 I < 2.74. 

(ii.) The scalar potential should not have an alternative minimum in which S gets a vacuum expectation value. A 
sufficient condition to ensure this is to set > 0 in Eq. Q. 

{in.) The electroweak oblique parameters S, T, and U [29], receive contributions from electroweak gauge boson self¬ 
energy diagrams involving the states in E. The measured oblique parameters constrain IA 3 I ~ 2 |A 4 | and favour 
A 3 < 0. This constraint causes the state to tend to cluster with the heavier scalars, as can be seen in the 
first and third plots in Fig. This will be the case unless the overall spectrum is highly compressed. 

{iv.) The SM Higgs boson decay width to two photons receives contributions from the charged scalars in E running in 
the loop. The LHG measurement of the h ^ 77 rate constrains A 2 more strongly than the unitarity constraint 
when 171(^0,r < 500 GeV. Since A 2 does not appear explicitly in the mass formulae in Eq. 0 or affect any of the 
electroweak production and decay rates, we satisfy this constraint by setting A 2 = 0 for tnis work. We checked 
numerically that this choice does not further restrict the range of the remaining parameters beyond the previous 
three constraints. 
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m^o,. [GeV] 

(n = 6 ) 

> II 

80 

76.98 

39.29 

120 

114.38 

55.75 

180 

170.48 

80.43 

300 

282.68 

129.81 


TABLE I: Maximum mass splitting between and — m^o,r, for particular values of in 

the two models. Values were obtained numerically by applying the theoretical and indirect experimental constraints 

described in Sec. 


rn^or = 80 GeV = 120 GeV = 180 GeV 





/n^o,r = 300 GeV 



FIG. 3: Parameter region allowed by the theoretical and indirect experimental constraints from Ref. [ 12 ] in the 
R 1 -R 2 plane, for four values of m(^o,r in the n = 6 model. The allowed region is shaded in yellow. Contours of 
constant A 3 are shown as solid curves with pink labels and contours of constant A 4 are shown as dashed horizontal 
lines with blue labels. Curves are labelled with their respective A 3 or A 4 value. 


In this study, it is most convenient to reparameterize the models in terms of physical masses and mass splittings. 
We take as inputs m^o,r as in Eq. Q, and replace A 3 and A 4 with the mass splitting parameters Ri and R 2 : defined 
as 


_ Amn ^ 171(^0,r ^ ^ 171(^0,i - 171(^0,r 

^ Amo m^o,i — m^o,r ’ ^ Am^^^ (m^^o,^ — m(^o,r)^^^ 

Here Ri parameterizes the mass of the highest-charged state, terms of the mass splitting Am^ = —m(^o,r 

as a fraction of the mass splitting between the lightest and heaviest state of E, Amo = rn(^o,i —m^o,r (see Fig.j^, while 
R 2 parameterizes Amo as a fraction of the maximum such splitting Am^^^ allowed after the theoretical and indirect 
experimental constraints are imposed on the model. The numerical values of Am^^^ allowed by the theoretical and 
indirect experimental constraints are given in Table [J for the four values of m^o,r used in our simulation. We study 
four mass slices, with m(^o,r = 80, 120 , 180, and 300 GeV. (LHC 8 does not provide any exclusions for m^o,r larger 
than 300 GeV in these models.) 

The mass splitting parameters are normalized to fall in the ranges Ri G [0, 1] and R 2 G [0, 1]. From Eq. Q, it is 
clear that as A 3 ^ + 2 |A 4 |, Ri -7- 0, and as A 3 ^ — 2 |A 4 |, Ri -7- 1. Furthermore, as illustrated in Fig. the mass 
splittings among the states are maximized when Ri ^ 0.5, while the spectrum collapses into two tightly-clustered 
groups of states when Ri -> 0 or 1 . Also, because R 2 parameterizes the overall mass splitting between the neutral 
real and imaginary scalars, the entire spectrum becomes compressed as R 2 0 , whereas the splitting between the 
heavier and lighter states is maximized when R 2 ^ 1 - 

The parameter regions allowed by the theoretical and indirect experimental constraints studied in Ref. [ 12 ] are 
plotted in Figs. 1^ and in the R 1 -R 2 plane, for the four m(^o,r values that we study. The yellow shaded regions are 
allowed. On eacn plot we also show the contours of constant A 3 and A 4 . 
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FIG. 4: As in Fig. but for the n = 8 model. 
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FIG. 5: Dominant decay modes of in the n = 6 model. Shown are ^ (orange/red) and ^ 

(purple/black). The decay ^ is severely kinematically suppressed over the entire parameter space and 

never dominates. The darker shade of each colour indicates that the gauge boson is emitted on-shell, while the 
lighter shade indicates it is off-shell. The region below the yellow curve is allowed by indirect and precision 

electroweak constraints. 


III. PREFERRED DECAY MODES 


The mass splittings and preferred decay modes between the physical scalars in our models vary substantially over 
the allowed parameter space. Understanding the collider constraints requires an understanding of the branching 
fractions of the unstable states over the R 1 -R 2 plane. In this section, we provide insight into the experimental results 
through examining the primary decay modes of the heavier branch of states over the full R 1 -R 2 plane, for m^o,r = 80, 
120, 180, and 300 GeV. Decays in the n = 6 model are shown in Figs. mzi and decays in the n = 8 model are shown 
in Figs. [SHTT] We do not plot the decay branching fractions of the lighter states, as these states always cascade down 
to through W emission as illustrated in Fig. 

Except for small regions in the transition from one dominant decay mode to another, the dominant mode has a 
branching fraction of at least 50%. We note that decays of a heavy state to another state in the heavier branch, such 
as ^ are never dominant. This is because they are highly suppressed by the lack of available phase space, 

making them effectively negligible over the entire parameter space. 


Our numerical results were obtained as follows. We implemented the two models in FeynRules version 2.0.26 m 
and generated the corresponding Universal FeynRules Output (UFO) files. We used these with MadGraph5_aMC@NL0 
version 2.1.2 [3T1[32] to compute the decay branching ratios of all the new scalars in a grid over the R 1 -R 2 plane for 
each of the four m^o,r slices. 

The branching ratios vary over the parameter space due to variation of the scalar mass splittings and the mixing 
angles ag, which in turn control the gauge couplings involved in the decays. The variation in the dominant decay 
mode with Ri is primarily due to variation in the decay couplings, which leads to the essentially vertical stripes in 
Figs.|5-fT The dependence on the mass splittings appears in regions of parameter space where one or more decays 
are transitioning between off-shell and on-shell. For example, in Fig. the ^ decay goes on-shell at a 
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FIG. 6: Dominant decay modes of H 2 in the n = 6 model. Shown are H 2 (yellow), H 2 Z 

(orange/red), and H 2 (purple/black). The decay H 2 H 2 ^W~ is severely kinematically suppressed 

over the entire parameter space and never dominates. The darker shade of each colour indicates that the gauge 
boson is emitted on-shell, while the lighter shade indicates it is off-shell. The region below the yellow curve is 

allowed by indirect and precision electroweak constraints. 



FIG. 7: Dominant decay modes of in the n = 6 model. Shown are (yellow), 

(orange/red), and (purple/black). The darker shade of each colour indicates that the gauge boson 

is emitted on-shell, while the lighter shade indicates it is off-shell. The region below the yellow curve is allowed by 

indirect and precision electroweak constraints. 


lower R 2 value than the other decay ^ leading the former to encroach into the parameter region that would 

have otherwise been dominated by the latter. This effect also shows up in Fig. because m^o,r < m^++, the decay 

of H 2 to goes on-shell while the decay to is still off-shell, leading to a small region of parameter 

space in which the former branching ratio surpasses the latter (the small island of black in the yellow region near the 
left side of the plot for m(^o,r = 120 GeV), despite the relative size of the couplings favoring decays to the latter 
To illustrate the interplay of coupling strengths and mass splittings in the pattern of branching ratios, we examine 
in more detail two slices through the parameter space of Figs. and In Fig. we show a slice through the 


= 180 GeV panel of Fig. |8| (decays of in the n = S model) at R 2 = 0.562. The relevant mass splittings 
normalized to the mass of the daughter gauge boson) are nearly constant with i?i, as are the 


12 


m^o 

(left plot in Fig 

and couplings (middle plot in Fig. 12). Instead, the variation in the branching ratios is due to 

the variation in the coupling, which passes through zero around Ri 0.6. This behavior can be traced 

back to the variation of the mixing angle a ^ . T his results in a big dip in the ^ branching ratio in favor 

of the decay ^ (right plot in Fig. 12). For small Ri where the decay to does dominate over 

the coupling is still smaller than the coupling; the larger branching fraction is instead due to the 


The appearance of these islands depends somewhat on the grid spacing of our scan. 
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FIG. 8: Dominant decay modes of in the n = 8 model. Shown are ^ (orange/red) and ^ 

(purple/black). The decay ^ is severely kinematically suppressed over the entire parameter space and 

never dominates. The darker shade of each colour indicates that the gauge boson is emitted on-shell, while the 
lighter shade indicates it is off-shell. The region below the yellow curve is allowed by indirect and precision 

electroweak constraints. 



FIG. 9: Dominant decay modes of H 2 in the n = 8 model. Shown are H 2 (yellow), H 2 Z 

(orange/red), and H 2 (purple/black). The decay H 2 H 2 ^W~ is severely kinematically suppressed 

over the entire parameter space and never dominates. The darker shade of each colour indicates that the gauge 
boson is emitted on-shell, while the lighter shade indicates it is off-shell. The region below the yellow curve is 

allowed by indirect and precision electroweak constraints. 


larger available phase space for the former mode. Note also that, although the coupling is significantly 

stronger than either of the others, the resulting decay is severely kinematically suppressed and never dominates. This 
is a generic feature of our models: due to the pattern of mixings among the isospin eigenstates, the gauge couplings 
between pairs of scalars tend to be largest when their mass splittings are smallest. The exception is the 
coupling, which does not depend on any mixing angles. 

we show a slice through the m^o,r = 120 GeV panel of Fig. M (decays of in the n = 6 model) at 
The three-band structure in Fig. [^emerges in a natural way from the interplay between couplings and 


In Fig. 


13 


R 2 = 0.429. 

mass splittings. At small Ri the coupling is sizable and the available phase space is large, so this decay 

dominates. As Ri increases, the mass of increases (see also Fig. [^, and the branching fraction 

drops due to the squeezing of the phase space, in spite of the growth of the coupling. Meanwhile, the 

coupling is passing through zero around Ri 0.3. This allows the decay to H^^Z to dominate for Ri 
between about 0.3 and 0.7. For high values of i?i, the rising coupling and its larger phase space allow 

the decay to to become dominant. 

Finally, we note that the maximum allowed mass splitting — m(^o,r is considerably smaller in the 

n = 8 model than in the n = 6 model (see Table |^. While on-shell decays first appear in the m^o,r = 120 GeV plots 
for the n = 6 model, they only appear in the m^o,r = 300 GeV plots for the n = 8 model. 
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FIG. 10: Dominant decay modes of in the n = 8 model. Shown are (yellow), Z 

(orange/red), and (purple/black). The decay H 2 ^W~ is severely kinematically suppressed 

over the entire parameter space and never dominates. The darker shade of each colour indicates that the gauge 
boson is emitted on-shell, while the lighter shade indicates it is off-shell. The region below the yellow curve is 

allowed by indirect and precision electroweak constraints. 



FIG. 11: Dominant decay modes of in the n = 8 model. Shown are (yellow), 

(orange/red), and H^^W~ (purple/black). The darker shade of each colour indicates that the gauge boson 

is emitted on-shell, while the lighter shade indicates it is off-shell. The region below the yellow curve is allowed by 

indirect and precision electroweak constraints. 


IV. COLLIDER STUDY 

LHC constraints on models with stable neutral particles (i.e., WIMP dark matter candidates) typically arise from 
searches for excess events with large missing transverse energy (^'t) recoiling against some combination of visible 
matter. These types of signatures have been studied both in the context of supersymmetry and in generic mono-jet 
searches by both ATLAS and CMS using LHC 8 data. Searches for jets+ ^ 7 ^ are typically only constraining for signals 
produced via QCD-strength interactions because of the large W/Z + jets backgrounds; nevertheless, we include them 
here for completeness. Signatures involving leptons + ^ 7 ^, for which the SM backgrounds are smaller, are typically 
more constraining for models in which the new physics is produced with electroweak-strength cross sections, as in our 
models. These leptons + signatures will be responsible for the direct-search constraints on our models. We take 
advantage of the large number of possible pairs of scalars that can be produced via 5 -channel 7 , Z, or W exchange in 
our models by simulating all such processes and combining the resulting signal events. 

The following searches were reproduced: 

• ATLAS opposite-sign dileptons with and no jets m-, 

• ATLAS trilepton plus [TTl [18] (the two studies employ distinct methodologies, and are reproduced sepa¬ 
rately); 

• ATLAS four or more leptons HU; 
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FIG. 12: Scalar mass splittings normalized to the emitted gauge boson mass (left), coupling strengths squared 
(middle), and branching ratios (right) for as a function of Ri in the n = S model for m(^o,r = 180 GeV and 
i ?2 = 0.562 (compare Fig. |^. Shown are the decays to (solid black curves), (dashed red curves), and 

(dotted blue curves). Note that the squared coupling (solid black curve) touches zero near 

Ri - 0.6. 






R^ 



FIG. 13: Scalar mass splittings normalized to the emitted gauge boson mass (left), coupling strengths squared 
(middle), and branching ratios (right) for as a function of Ri in the n = 6 model for m^o,r = 120 GeV and 

i ?2 = 0.429 (compare Fig. [^. Shown are the decays to (dashed red curves), Z (solid black curves), and 

(dotted blue curves). Note that the H^W~ squared coupling (dotted blue curve) touches zero near 

Ri - 0.3. 


• ATLAS dileptons with razor variables [20] : 

• ATLAS hadronic di-r plus [22]; 

• ATLAS same-sign dileptons plus jets ED; 

• ATLAS monojet [2^1^ : 

• ATLAS multi-jets plus [25] : 

• CMS multi-leptons (dilepton, trilepton, multi-lepton) with [26] : 

• CMS monojet [28] . 

We performed the simulation following the methodology of Ref. [33]. For each of the four mass slices, m^o,r = 80, 
120 , 180, and 300 GeV in each model, we generated parameter points on a grid in Ri and R 2 , focusing on the region 
that is allowed by the theoretical and indirect experimental constraints. We also simulated a few additional points 
for added sensitivity in regions where the experimental sensitivity was varying rapidly. The parameter points that we 
simulated are shown as green dots on our exclusion plots. 

We do not simulate parameter points for which Amo = — m^o,r < 10 GeV. As we will see, this compressed- 

spectrum region is particularly difficult to detect because the leptons resulting from the decays tend to be too soft to 
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pass the selection cuts. Indeed, we will find that parameter regions with Amo ^ 20 GeV cannot be excluded using the 
LHC 8 searches that we recast. This parameter region could be picked up by jets + searches, but we find that the 
scalar pair production cross sections are too small for LHC 8 to have any sensitivity in this channel. Very small mass 
splittings could result in macroscopic decay lengths for charged scalars; while dedicated searches for such signatures 
exist, we will not consider them here. 

Using MadGraph5_aMC@NL0 and our FeynRules model files, we generated 50,000 events for each of the possible pairs 
of scalars that can be produced via 5-channel and 7 */^ exchange, in association with zero or one hard jet(s). 
This amounted to 32 possible scalar pairings in the n = 6 model and 44 in the n = S model. Each simulated event 
sample was then matched and merged, decayed, showered, and hadronized using Pythia 6 [34]. The events were then 
passed through the Delphes 3 [35| detector simulation multiple times: for each of the experimental searches that we 
recast, we adjusted the identification efficiencies and jet algorithm settings in Delphes 3 to match the working point 
used in the experimental analysis. We then used the Seer [36| analysis program to apply the trigger and kinematic 
cuts of each experimental analysis and determine the total signal cross section from all production processes that fall 
into each of the signal regions. Additional Gaussian smearing of the beyond that already present in Delphes 3 
was introduced via Seer to account for the effect of pileup m This was necessary in order to reproduce the cut-flow 
tables for each of the experimental searches. 

At each simulated parameter point, we select as most sensitive the experimental search that yields the largest value 
of logio(A^sig/V 95 ), where Vgig is the number of signal events computed from the cross section in the signal region 
for the appropriate integrated luminosity and V 95 is the corresponding 95% confidence level exclusion threshold from 
the experimental analysis. If log;Lo(^sig/A' 95 ) > 0 for the most sensitive search, the parameter point is excluded. 
We do not combine significances from different search channels—to do so would require knowledge of the statistical 
correlations among the various searches which we do not have. In between parameter points, we interpolate linearly 
in the value of log]^o(^sig/A^ 95 ) found for the most sensitive search at each point. 

We find that the most sensitive analyses are the ATLAS dilepton, trilepton and four-lepton searches. Several other 
searches also exclude parameter space for low values of m^o,r, but this is primarily due to the large pair production 
rates when the scalar masses are low, which overcome the inherently smaller acceptances and/or larger backgrounds 
in these searches. 

Our results are shown in Fig. [^for the n = 6 model and Fig. [^for the n = 8 model. For each of the mass slices 
m^o,r = 80, 120 , 180, and 300 GeV, we show the excluded region of the R 1 -R 2 plane (darker-colored regions above 
and to the left of the solid black line). The dashed black line indicates how the excluded region would expand if the 
signal cross sections were all increased by a /c-factor of 1 . 2 . As before, the parameter region below the thick teal curve 
is allowed by the theoretical and indirect experimental constraints. In each plot the thick horizontal light pink line 
indicates the R 2 value at which Amo = — m(^o,r = 20 GeV, below which the spectrum is highly compressed. 

We now describe the results for each model. 


A. Results for the n = 6 model 


In the n = 6 model, the parameter space is strongly excluded at low m(^o,r so long as Amo > 20 GeV. As m(^o,r 
increases, the excluded region shrinks until only a small sliver of parameter space around R 2 ^ 0.3 is excluded for 
171(^0,r = 180 GeV. LHC 8 makes no exclusion at m^o,r = 300 GeV or above. 

For 171(^0,r = 80 GeV, shown in the upper left panel of Fig. 
cuts ranges from ^ 30 pb for small R 2 to ^ 16 pb for large R 2 . 
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the total scalar pair production cross section before 
The ATLAS trilepton+^ 7 ^ searches [T7l[T8] provide the 


most sensitive exclusion for all parameter points scanned. Due to the large signal cross sections, even small acceptance 
rates result in large numbers of events in the experimental signal regions, resulting in the strong exclusion as indicated 
by dark colors. The excluded region ends abruptly at Amo — 20 GeV. For R 2 values below this boundary, the highest 
lepton pt is rarely above 10 GeV, unless the scalars are produced boosted or in association with an energetic jet, 
which is rare. The ATLAS and GMS offline leptonic triggers are typically set at pr > 10 GeV, resulting in a drastic 
plunge in acceptance. Similarly, the tends to be low when the spectrum is compressed; most (but not all) of 
the searches require > 50 GeV, further reducing sensitivity to the compressed-spectrum region. This kinematic 
acceptance boundary is apparent in many of the exclusion plots for both the n = 6 and n = 8 models. 

For m^o,r = 120 GeV, shown in the upper right panel of Fig. 
cuts ranges from ^ 8.2 pb for small R 2 to 
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the total scalar pair production cross section before 
3.5 pb for large R 2 . The ATLAS trilepton + searches [T7l[T8] provide 


the most sensitive exclusion for R 2 ^0-5, in particular the search region with same-flavor opposite-sign dilepton mass 
smaller than the Z mass, which picks up the off-shell Z boson decays. For R 2 > 0.5 and Ri > 0.75, the most sensitive 
search is the ATLAS four-lepton search m, due to the smaller backgrounds and the predominance of ^ 


and H 2 


Z decays in this region (see Figs. [^and[^. The unexcluded region in the upper right corner of this 
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FIG. 14: Combined exclusion regions in the n = 6 model from all examined LHC8 analyses. The darker-colored 
regions above and to the left of the solid black line are excluded at 95% confidence level by at least one analysis. 
The dashed black line indicates how the excluded region would expand if the signal cross sections were all increased 
by a /c-factor of 1.2. The parameter region below the thick teal curve is allowed by the theoretical and indirect 
experimental constraints. The thick horizontal light pink line indicates the R 2 value at which 

Amo = Tn(^o,i — 171(^0,r = 20 GeV. 


plot results from the loss of the four-lepton signatures as the decays involving W emission go on shell (the resulting 
two-lepton final states suffer from larger backgrounds). 

For m^o,r = 180 GeV, shown in the lower left panel of Fig. 14, the total scalar pair production cross section before 
cuts ranges from ^ 1.7 pb for small i ?2 to ^ 0.7 pb for large R 2 . The ATLAS trilepton + searches [ITl [18] are 
responsible for the exclusion between R 2 values of about 0.2 and 0.3, while the ATLAS four-lepton search m is 
responsible for the thin strip of excluded parameter space near Ri ^ 0.75. Within the parameter region allowed by 
theoretical and indirect experimental constraints, only a small sliver of parameter space around R 2 ^ 0.3 is excluded 
by the ATLAS trilepton + searches. The largest contribution to the trilepton + signal is from the 
pair production mode, which decays in this region of parameter space mainly to Z, followed by very soft 

decays of the states. As R 2 increases from 0.3 to 0.4, the decay transitions to due to the nearby 

kinematic threshold (see Fig. so that the third lepton is lost. The production cross sections for the heavier scalars 
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FIG. 15: Combined exclusion regions in the n = 8 model from all examined LHC8 analyses. The darker-colored 
regions above the solid black line are excluded at 95% confidence level by at least one analysis. The dashed black 
line indicates how the excluded region would expand if the signal cross sections were all increased by a /c-factor of 
1.2. The parameter region below the thick teal curve is allowed by the theoretical and indirect experimental 
constraints. The thick horizontal light pink line indicates the R 2 value at which Atuq = m^o,i — m(^o,r = 20 GeV. 
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the total scalar pair production cross section 


also decrease with increasing R 2 . 

For 171(^0,r = 300 GeV, shown in the lower right panel of Fig. 
before cuts ranges from ^ 0.2 pb for small R 2 to ^ 0.08 pb for large R 2 . After accounting for branching ratios and 
acceptances, it is clear that there is simply insufficient luminosity collected at LHG8 to result in any sensitivity to the 
model for this and higher masses. 


B. Results for the n = 8 model 

In the n = S model, the parameter space at low m(^o,r is again excluded so long as Amo ^ 20 GeV. The excluded 
region shrinks rapidly with increasing m^o,r until LHG8 can make no exclusion at m^o,r = 180 GeV or above. 

The most important difference between the two models is that the maximum mass splitting between the lightest 
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and heaviest scalars, is considerably smaller in the n = 8 model than in the n = 6 model, 

by about a factor of two (see Table ]^. Therefore, a much larger fraction of the allowed n = 8 model parameter space 
lies in the compressed region with Amo < 20 GeV (below the thick light pink horizontal line in Fig. 15). The smaller 
mass splittings over the entire n = 8 parameter space result in softer leptons, jets, and hence overall lower 

acceptance rates in the LHC 8 analyses. This results in weaker exclusions in the n = 8 model in spite of the larger 
number of scalars, their higher weak and electric charges (leading to more like-sign leptons from longer decay chains), 
and the correspondingly higher total scalar pair production cross sections before cuts. The emitted leptons are just 
too soft. 

For 171(^0,r = 80 GeV, shown in the upper left panel of Fig. the total scalar pair production cross section before 
cuts ranges from ^ 90 pb for small R 2 to ^ 62 pb for large i? 2 ^rhe entire parameter space with R 2 ^ 0.5 is excluded. 
Similarly to the n = 6 model, the excluded region ends around Amo — 20 GeV, below which the highest lepton pt is 
rarely above 10 GeV, so that most events fail the offline leptonic triggers. The most sensitive exclusion comes from 
the published ATLAS trilepton + study [18] . This study captures a broader range of signal characteristics than 

the corresponding conference note and includes signal regions that are more sensitive to low dilepton invariant 
masses (called the SR^ra bins 1-4 in Ref. [E]). In particular, the best constraints on the n = 8 model come from 
SROra bin 2 , in which the invariant mass of the same-flavor opposite-sign lepton pair lies between 12 and 40 GeV, 
the transverse mass of the remaining lepton and the lies below 80 GeV, and > 90 GeV. The most constraining 
bin in Ref. m is similar, but with a softer cut on ^ 7 ^, which results in larger background. The ATLAS four-lepton 
search m also redundantly excludes part of the parameter space, but loses sensitivity for larger values of Ri, where 
the dominant decay of changes from Hf^Z to H^W^. 

For m^o,r = 120 GeV, shown in the upper right panel of Fig.f^ the total scalar pair production cross section before 
cuts ranges from ^ 20 pb for small R 2 to ^ 12 pb for large R 2 . The excluded region is understood similarly to the 
previous mass slice, with the same mechanisms at work. When the smaller cross sections are combined with branching 
ratios and acceptance rates, the signal events in the published ATLAS trilepton + selection [T8| wind up falling 
into multiple signal bins, rather than being as concentrated in a single bin as in the m(^o,r = 80 GeV slice. Since we 
do not have information on the correlations of the systematic uncertainties among different bins, we cannot employ 
the full usefulness of the CLq method to combine the bins to extend the exclusion boundaries. In addition, decays to 
pairs of off-shell Z bosons are common, which reduces the trilepton signal from WZ production. A four-lepton search 
is intrinsically less sensitive than a trilepton search due to the smaller leptonic branching ratio of the Z compared to 
the W. 

For m^o,r = 180 GeV, shown in the lower left panel of Fig. the total scalar pair production cross section before 

cuts ranges from ^ 4.7 pb for small R 2 to ^ 2.7 pb for large None of the scanned parameter space is excluded by 

LHG 8 searches. A small region of parameter space at R 2 > 0.9 and Ri ^ 0.7 would be excluded if a signal /c-factor 
of 1.2 were included. The sensitivity in this region comes from the ATLAS opposite-sign dileptons + search with 
no jets m, which is sensitive to the off-shell Z decays from the and that are produced in association with 
Hi states. Jets from the decays of the Hi states are soft enough that they evade the jet veto in the search. 

For 171(^0,r = 300 GeV, shown in the lower right panel of Fig. the signal rate after cuts in any channel is less 
than 10 % of that required for an exclusion over the entire scanneo^arameter space. 


V. DISCUSSION AND CONCLUSIONS 

We studied the direct-search constraints from LHC8 data on models in which the Higgs sector is extended by a Z 2 - 
odd complex scalar electroweak multiplet S with isospin T = 5/2 (n = 6) or 7/2 (n = 8) and the same hypercharge as 
the SM Higgs doublet. These models can be probed by recasting the dedicated searches for supersymmetric particles 
performed by the CMS and ATLAS collaborations. We showed that, even after imposing theoretical and indirect 
experimental constraints (mainly from the electroweak oblique parameters), the LHC8 searches provide nontrivial 
further constraints on the parameter space of the two models. 

The data from LHC8 excludes the majority of the remaining parameter space for m^o,r = 80-120 GeV in both the 
n = 6 and n = 8 models, except for parameter regions in which the entire spectrum is compressed within a mass 
splitting Amo < 20 GeV. The loss of sensitivity in the compressed spectrum region is caused by the extreme softness 
of the cascade decay products, leading most of the signal events to fail the trigger and minimum energy requirements 
for lepton tagging. The LHC8 constraints rapidly disappear for m^o,r > 120 GeV, until there is little sensitivity above 
m^o,r = 180 GeV. These results are in general agreement with the LHC8 constraints on the electroweak gauginos 
in supersymmetric models studied in Ref. [33| , in which the same collection of LHC8 searches showed sensitivity to 
masses of the lightest supersymmetric particle primarily smaller than 100 GeV. 

The LHC8 constraints on our models are much weaker than those on the scalar septet extension of the SM Higgs 
sector [38| studied in Ref. [39]. Using similar searches from ATLAS and CMS, Ref. [39] found that septet masses 
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below ^ 400 GeV were robustly excluded, and that multi-lepton searches provided the strongest constraints. The 
dramatically enhanced LHC 8 sensitivity to the septet model compared to our sextet and octet models is due to the 
nonzero vacuum expectation value carried by the septet and the mixing between the neutral scalar in the septet model 
and the SM Higgs boson in that model. These features allow the lightest septet state to decay into SM particles, 
particularly into energetic gauge boson pairs, leading to final states with two to four high-energy leptons. The mass 
reach of the LHC 8 exclusion is then controlled by the decrease in the production cross section with increasing scalar 
masses. 

With Run 2 of the LHC underway at 13 TeV center-of-mass energy, the next round of searches will likely provide 
stronger constraints on our models, extending to higher scalar masses. The pair production cross section for fixed 
scalar masses grows with increasing proton-proton center-of-mass energy due to the rise in the antiquark parton 
density. Higher integrated luminosity will also aid the searches, but with diminishing returns. The analyses that 
we use compare an observed event number to theoretical predictions for the background cross section; unless these 
backgrounds can be normalized from data, such searches suffer from systematic uncertainties that ultimately limit 
the benefit of greater luminosity. In addition, higher instantaneous luminosity will necessitate an increase in trigger 
thresholds, which will further reduce the sensitivity to the compressed spectrum regions. Discovery prospects at a 
higher energy proton-proton collider would face similar issues. 

Finally, we would be remiss not to mention the fact that the scalars in S could be pair produced via 5 -channel 
photon and Z exchange at a lepton collider such as the International Linear Collider (ILC) [40]. Such a machine 
would be able to probe scalar masses up to the pair production threshold ~ \/ 5/2 with small backgrounds. ILC search 
prospects in the compressed spectrum region would require a dedicated study. 
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Appendix A: Masses and mixing angles 


In this section we give some of the mathematical details used in the derivation of the mass spectrum and mixing 
angles in Sec. [ll| 


Eor a complex scalar multiplet E with hypercharge T = I (normalized so that Q = + T/2), the most general 
gauge-invariant and Z 2 -invariant renormalizable scalar potential was given in Eq. Q, in which 4> = and 

E = CE* are the conjugate multiplets. Here is the second Pauli matrix and the conjugation matrix C for the large 
multiplet is an anti-diagonal n x n matrix. Eor n = 6 and 8 , the matrix C is given by 


^(n=6) 


/ 0 0 0 0 0 I\ 

0 0 0 0 -1 0 

0 0 0 1 0 0 

00-1000 

0 1 0 0 0 0 

V-I 0 0 0 0 0/ 


C(n=8) 


/ 0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 1 
0 0-10 
0 10 0 
V-I 000 


0 0 0 I\ 
0 0-10 
0 10 0 
-10 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 / 


(AI) 


Taking A 4 real and working in unitarity gauge, the term involving A 4 in 


the scalar potential of Eq. § reduces to 


A4$V“$I]tT“E + h.c. 


-Xiih + vf 


E^T“E + E^T+E 


(A 2 ) 
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where ± zT^. For n = 6 the generators T“ are given by 


T"*" = 

(n=6) 


/O v/5 0 0 0 0 \ 

0 0 2 V 2 0 0 0 

0 0 0 3 0 0 

0 0 0 0 2 V 2 0 

0 0 0 0 0 

\o 0 0 0 0 0 y 


(^(n=6)) 


while for n = 8 they are 


T+ 

^ (n=8) 


rfn= 6 ) = ^diag(5, 3, 1, -1, -3, -5), 


/O \/7 0 0 0 0 

0 0 2\/3 0 0 0 

0 0 0 \/l5 0 0 


0 0 
0 0 
0 0 
0 0 
Vo 0 


0 0 \ 

0 0 

0 0 

0 4 0 0 0 

0 0 yi5 0 0 

0 0 0 2 v^ 0 

0 0 0 0 y? 

0 0 0 0 0 / 


(^(n=8)) ’ 


rfn=8) = 2 diag (7, 5, 3, 1, -1, -3, -5, -7). 


(A3) 


(A4) 


(A5) 


(A 6 ) 


The terms T^T S and S’f’T+S split the masses of and and cause mixing between states with the same 
electric charge but different isospin. For n = 6 or 8 , these two terms can be written as 

n/2 

I]tT“E = -(-l)"/2+yo*^o* + A/n2 - 4Q2 (_i)"/2+Q+i 
2 Q =1 

n/2 

S’l'T+S = -(-l)”/2+yo^o + ^ - 4(32(_i)"/2+Q+i ,^+<3,^-0. (A7) 

^ Q =1 


Writing the neutral state in terms of its real and imaginary components, = (C°’^ + *C°’*)/a/ 2, we find a mass 
splitting between the components, 


m^o.r = + -v^ 

m^o.i = M"^ + 


A 2 + yA3 + ;(-l)”/2+iA4 


A 2 + iA3 + |(-l)”/% 


= M 2 + -z; 2 A„, 


rri/- 




(A 8 ) 


The mass matrices for the pairs of scalars with positive electric charge Q = 1, • • •, f — 1 are given in the basis 

by 


_ / M2 + lu2(4A2-(2g-l)A3) iu2A4vA)^^4^(-l)”/2+«+iV 

- 4^ (-l)”/2+Q+i M2 + |t;2(4A2 + (2Q + l)A3) y ’ ^ ^ 

which we diagonalize to find the mass eigenvalues, 

(^^2 + ^Aa T ^ ]/ Q^A| + (n2 - 4Q2 )a2^ 

= TO^o.r + ^n(-l)"/^A4 =F ^Q^A| + (n2 - 4g2)A|^ . (AlO) 
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The mass eigenstates and are defined in terms of the weak eigenstates by Eq. ^ such that is the lighter 
state and is the heavier state. The mixing angle aq G [—f, f] is given by 


“ ' v/n» - 402 A 4 




\Jrfi - 4Q^ A 4 


QA 3 + v/Q2^i + {n2-4Q2)A2‘ 
There is only one state with Q = nl2. Its mass is given by 


m ^„/2 = 


+ ( 4 A 2 - (2Q - 1)A3) = m^o.r - (^A3 + 2(-l)"/^+^A4j . 


(All) 


(A12) 


Appendix B: Details of the re-parameterization 

For this study, it is convenient to describe the parameter space in terms of the mass differences Amo = '^(^o,r 
and Amn = since these two mass splittings are monotonic in —A 3 and IA 4 I. In terms of the original 

Lagrangian parameters, the mass splittings are given by 


/ Ti 

Amo = — m(^o,r = Wm^o,r + —(—l)5'r’^A4 — m(^o,r , 

/ Tl 

mQ0,r = Jm^o.r- “ [A 3 + 2 (-l)t+iA 4 ] - m^o,, 


Am- = TO^+n 

We defined normalized versions of these splittings, each lying in the range [0, 1], as 

Amo m^o,i —171(^0,, 


_ Amn m^+n- 171(^0,r 
Amo m^o,i — m^o,r 


R2 = 


Am^^^ (m(^o,i — m^o,r)^^^ 


(Bl) 

(B2) 

(B3) 


The maximum value of Amo allowed by the theoretical and indirect experimental constraints (mainly from the 
electroweak oblique parameters), Am^^^, can be written parametrically for the n = 6 and n = S models as 


Amn 


2.2 T 0.94m(^o,r 

mt;o,r < 530 GeV, 

(B4) 

809.4 — 0.72171(^0,r + 0.00024m^o,r 

m^o.r- > 530 GeV, 

6.37 + 0.4Im(^o,r 

niQo.r < 809 GeV, 

(B5) 

687.7 — 0.57m^o,r + 0.000I7m^o,r 

> 809 GeV. 


The numerical values of Am^^^ for our chosen mass slices were given in Table 


The relations in Eqs. (B3), (BI), and (B 2 ) can be inverted to obtain A 3 and A 4 as follows: 


A3 = 


4i?2 Am^^^ 


nv^ 


[2 (I - 2Ri) m^o,r + (I - 2Rl) R 2 Am^^^] , 


+ii, . 


(B 6 ) 

(B7) 


The physical masses can also be expressed as 

m^o,i = m^o,r + i?2 Am^"^^, 


I 


m^+Q = y ^^o,r + -i?2(Am^^^)2 + R2 
m^+^ = m^o,r + i?ii? 2 Am^'^^ , 


TTL(^0,r 


2n 


Am^^^, 


(B 8 ) 
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State 

Mass [GeV] 

Arrii^i-i [GeV] 


80.0000 

- 


80.1960 

0.196 


80.8320 

0.636 


82.1524 

1.3204 


89.2535 

7.1011 


89.8254 

0.5719 


90.0000 

0.1746 


TABLE II: Mass spectrum and splittings for a sample point in the compressed-spectrum region, Amo = 10 GeV. 
For this point, n = 6, m(^o,r = 80 GeV, Ri = 0.215 and R 2 = 0.130. Compare the charged pion mass at 0.140 GeV. 


where 


M= yj (n2 - 4Q2) (2m^o.r + + 4Q2 [2 (1 - 2Ri) + (1 - 2Rl) . 


Appendix C: Compressed spectra versus Pythia 

Our models include parameter regions in which the spectrum is highly compressed. In such regions, simulation 
problems can arise because the automatic parton-level computation of partial widths and branching ratios using 
MadWidth [32] does not take into account hadronization effects arising from decays involving a highly off-shell W or 
Z boson. 

For example, consider the case of decay for a parameter point at which the mass splitting 

m^+ — 171 (^ 0 ,r is larger than the mass of a tt^ but smaller than the mass difference of the tt^ — tt^ system, as in 

Table [nj Physically, a hadronically-decaying off-shell W will emerge in this case as a single tt^. However, MadWidth 
treats all decay products as bare partons and will decay the off-shell IT to a quark-antiquark pair. When Pythia [34] 
hadronizes the resulting quark-antiquark pair, it tries to create two pions, finds that this is forbidden by the available 
phase space, and discards the event. 

This is a problem because it removes generated events that could otherwise potentially be detected via, e.g., a 
hard initial-state jet radiation. We solve this problem by manually removing any hadronic IT branching ratio from 
the MadGraph param_card files in which the IT invariant mass is less than the mass of two correspondingly-flavored 
mesons and rescaling the remaining IT decay branching ratios to sum to one. 

Should we worry about changing the potential detector signature of the very soft IT boson? A simple kinematical 
argument indicates no. In the rest frame of the decaying scalar, the momentum of the outgoing decay products must 
satisfy, e.g., \P(^o,r \ + |pjets| ^ ~ above example. In this frame the decay product(s) of the IT will 

have momentum less than the scalar mass splitting. To make the IT decay product(s) visible in the detector, they 
must be boosted to an ultrarelativistic velocity; for example, a 10 GeV charged pion requires a relativistic boost factor 
7 70. This requires that the parent scalar itself be boosted by the same relativistic boost factor: for a scalar mass 

around 80 GeV, this implies a transverse momentum of more than 5 TeV, which is not even kinematically possible at 
the 8 TeV LHC. 
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